Quantum entanglement of nanocantilevers 
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We propose a scheme to entangle two mechanical nanocantilevers through indirect interactions 
mediated by a gas of ultra cold atoms. We envisage a system of nanocantilevers magnetically 
coupled to a Bose-Einstein condensate of atoms and focus on studying the dark states of the system. 
These dark states are entangled states of the two nanocantilevers, with no coupling to the atomic 
condensate. In the absence of dissipation, the degree of entanglement is found to oscillate with time, 
while if dissipation is included, the system is found to relax to a statistical mixture of dark states 
which remains time independent until the inevitable thermal dephasing destroys the nanocantilever 
coherence. This opens up the possibility of achieving long-lived entangled nanocantilever states. 

PACS numbers: 42.50.Dv, 37.30.+i, 03.65.Ud, 85.85.+j 



I. INTRODUCTION 

There have been long-standing arguments regarding 
the validity of quantum mechanics in the macroscopic 
world. The last two decades have seen an unprecedented 
rise of interest in studying the crossover between quan- 
tum and classical mechanics. These studies, both theo- 
retical and experimental, have the potential to improve 
our understanding of decoherence, which is assumed to 
degrade any quantum system to its classical counterpart. 
Over the past few years, there have been many fascinat- 
ing experiments to prepare and detect quantum states of 
macroscopic objects [l|. In this quest for studying the 
level of "quantumness" present in macroscopic objects, 
tremendous progress has been achieved in exploring the 
uantum regime of micro- and nano-mechanical systems 
. These offer a very promising playground for studying 
the quantum-classical crossover. Miniaturized cantilevers 
contain a macroscopic number of atoms and can be fabri- 
cated to have very high resonant frequencies and exceed- 
ingly large quality factors [1], thereby guarding against 
the effects of decoherence. 

The advancement in techniques such as laser cooling of 
mechanical resonators has brought quantum state 

preparation within experimental reach. Experimental 
advances have also fueled a surge of interest in coupling 
these condensed-matter systems to other quantum sys- 
tems with well-understood quantum properties, such as 
Bose-Einstein condensates or electron spins. Such hy- 
brid quantum systems are promising candidates for ultra- 
precise measurements Very recently, O'Connell et 
al. 0] were able to cryogenically cool a mechanical res- 
onator to its quantum ground state, and were also suc- 
cessful in strongly coupling it to a superconducting qubit 
to read out the motion of the resonator. This success 
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heralds a new era in investigating the quantum behavior 
of nanomechanical systems. 

The study of entanglement and superpositions of 
macroscopic objects is of prime interest. These phenom- 
ena have been experimentally demonstrated for meso- 
scopic systems, by interferometry of molecules and 
entangling of atomic ensembles [llj. There are propos- 
als to generate spatially separated superposition states 
of nanomechanical systems using controlled interaction 
through a Cooper-pair box [l2l | , and to use nanomechan- 
ical systems as a quantum bus for quantum informa- 
tion processing [l^, [3- Bose and Agarwal have pro- 
posed a scheme to entangle two nanocantilevers through 
a Cooper-pair box [Tsj . while Treutlein et al. have pre- 
sented a scheme to entangle the vibrational mode of a 
nanocantilever with the spin degree of freedom of an ul- 
tracold Bose-Einstein condensate (BEC) [31 . Here, a 
nanocantilever is coupled to a cloud of ultracold atoms 
through an external magnet. This scheme is interesting 
since it would allow us to couple two nanocantilevers to 
the same cloud of ultracold atoms, through which the 
nanocantilevers can interact with each other. Hunger et 
al. were recently able to resonantly couple the vibrations 
of a micromechanical oscillator to the motional degree of 
freedom of an ultracold BEC in a trap using the atom- 
surface interactions In their setup, the ultracold 
BEC was trapped at a few /im from the surface of the 
cantilever. When the cantilever vibrates, the surface po- 
tential becomes time dependent and modulates the trap- 
ping potential for the ultracold BEC. This in turn excites 
the atomic motion and generates coupling between the 
two. 

Motivated by recent experimental results, we here the- 
oretically investigate the possibility of entangling two 
nanocantilevers coupled via an ultracold Bose gas. We 
model the nanocantilevers as quantum harmonic oscilla- 
tors which have been cooled near to their ground states, 
so that the average number of excitations is much smaller 
than unity. The ultracold Bose gas is assumed to com- 
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prise N two-level atoms, and is modeled using the Dicke 
formalism [l^. We study the time evolution with and 
without dissipation, and find that entanglement arises 
in either case. Unitary evolution of the system leads 
to generation of time-varying entanglement between the 
two nanocantilevers. The dynamics becomes increasingly 
complex the more highly excited the studied system is. 
In the presence of dissipation in the atomic gas, the two 
nanocantilevers will be entangled in a steady state, which 
is decoupled from the atomic gas and has a life time that 
is only limited by the inevitable dephasing caused by the 
environment. This opens the possibility of achieving long 
lived "macroscopic" entangled states. 

The paper is organized in the following way. Section HIl 
introduces the theoretical model and the physical setup 
for a system of two nanocantilevers interacting via an ul- 
tracold Bose gas. This forms the framework for Sec. IIIII 
in which the unitary evolution of the system is studied 
using the Schrodingcr and Heiscnberg time evolution in 
different excitation manifolds. Decoherence effects are 
included in the description in section ITVl while the entan- 
glement present in the system is quantified in section [V] 
Finally, a discussion concludes the paper in section IVTl 



II. THEORETICAL MODEL AND PHYSICAL 
SETUP 



We propose a scheme to generate entangled states of 
two nanocantilevers through indirect interactions medi- 
ated via an ultracold Bose gas. Let us assume that two 
identical nanocantilevers, each mechanically vibrating in 
its fundamental flexural mode [l^, have been precooled 
near to their ground states. We model the two nanocan- 
tilevers as quantum harmonic oscillators, with a^(a) and 
b^{b) as the creation (annihilation) operators for vibron 
excitations in cantilever a and 6, respectively. The ultra- 
cold gas is described as a collection of N two-level atoms, 
modeled by the Dicke formalism [3 with J+ ( J_) as the 
collective spin raising (lowering) operator. 

The physical setup is shown in fig. [TJ The two iden- 
tical nanocantilevers, separated by a distance 2d along 
the y-axis, are assumed to be fabricated on an atom chip 
with strong ferromagnets attached to their tips. Equidis- 
tant from the tips, at a distance d, an ultracold Bose gas 
is confined in a microtrap. The magnetic moment /x of 
each ferromagnet is pointing in the cc-direction and the 
ultracold atoms are subject to a stationary magnetic field 
of strength Bq along the z-direction, which becomes the 
quantization axis for the ultracold atoms. For the pur- 
pose of providing an identical trapping potential for all 
the hyperfine states of the atoms [la] we envisage an op- 
tical or clectrodynamic trap [2^. In such a trap all the 
atoms in the BEC couple simultaneously to the quantized 
motion of the resonator. 

Under the dipolc approximation, the x-component of 
the magnetic field at the center of the trap, produced by 



the ferromagnet on the tip of cantilever a, is 



Br 



4Tr[d + ya.{t)]^ 



(1) 



where fiQ is the vacuum permeability and ?/a(0 is the 
time dependent defiection of the tip of nanocantilever a. 
For small displacements of the cantilever this expression 
can be expanded to linear order in ya(t), so that 



Br,: = - 
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Vait) 



(2) 



Thus yait) transduces the vibrational motion of can- 
tilever a to an oscillating magnetic field given by 



B = Gmyait)x 



(3) 



at the location of the ultracold atoms, where 
G'm=3^/xo/47rc?^ is the magnitude of the magnetic field 
gradient in the y-direction. By making use of, e.g., the 
quadratic Zeeman effect selected hyperfine states of the 
atoms in the ultracold gas can be decoupled from other 
states. Transitions from \F = l,m = 0) = |0) to 
|F = l,m = —1) = |1) can, for instance, be resonantly 
coupled to the quantized bending motion of the cantilever 
by tuning the Larmor frequency. In this way coupling 
is achieved between the fundamental bending mode of 
cantilever a and the collective spin of the ultracold gas 
through the Zeeman interaction. 



B(t) = hKa 



1 



(J+ + J_)(a + at), (4) 



where fi atom is the magnetic moment of the ultracold gas 
and Ka=fJ-BGmV^ao/\/8h is the coupling constant for 
the collective atomic spin coupled to a single vibrational 
mode of cantilever a. Here TV is the number of atoms 
in the gas and oq = ^Jhllm^giOQ is the amplitude of the 
zero-point fluctuations of cantilever a, which has effective 
mass meg and angular oscillation frequency wo- The 
coupling between cantilever h and the ultracold Bose gas 
is described by an analogous term. 

In the regime of strong coupling between the collec- 
tive spin of the ultracold gas and the vibrational mo- 
tion of the nanocantilevers, the system is analogous to 
a conventional cavity-QED system. For a nanocan- 
tilever with mass ra ~ 10~^^ kg and resonance frequency 
ujQ /2tt ~1 MHz, which is separated from an atomic cloud 
of ~ 100 atoms by a distance of d ~ 250 nm, the in- 
teraction produced by a small [2l| disk-shaped magnet 
containing N^ag ^ 10® nickel atoms corresponds to a 
coupling constant Ka ^ 100 Hz. 

In the particular geometry we envisage for physically 
entangling the two nanomechanical systems, direct inter- 
action between the magnetic dipoles on the two cantilever 
tips can be neglected. To see this, one can compare the 
two interaction strengths, one due to the direct interac- 
tion between the magnetic cantilever tips and the other 
due to the interaction between either cantilever magnet 
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FIG. 1: (Color online) Physical setup for the proposed scheme 
for entangling two nanocantilevers. Two identical nanocan- 
tilevers have strong ferromagnets attached to their tips. The 
cantilevers are placed equidistant from an ultracold gas of 
atoms, which is confined to a microtrap. Each ferromagnet 
couples the vibrational motion of a nanocantilever to the col- 
lective spin J of the ultracold gas. 



tem is 

— = uja{J.+a^a + Pb) + [K-^{a + b)J+ + h.c.l6) 
ri \/ N 

= Ho + Hi, 

where Hq = u!o{Jz + a^a + b^b) and Hj = K{a + b)j+ + 
h.c. are the free and interaction parts of the Hamiltonian 
respectively, and is the z component of the collective 
angular momentum of the ultracold gas. In the Dicke 
model, the collective spin of the ultracold BEC is J = 
N/2, where N is the number of atoms. A general Dicke 
state is represented as \M, J), where J+M is the number 
of excitations and J — M is the number of atoms in the 
ground state [23| . In this notation, | — J, J) represents the 
ground state and \J, J) the highest excited state. Spin 
raising and lowering operators act on a general Dicke 
state as 

J+\M,J) = VJ{J + 1) - M{M + 1)\M + 1, J) (7) 
J-\M,J) = y/J{J + 1) - M{M - 1)|A/ - 1, J). (8) 

After these initial considerations, we proceed to study 
the unitary evolution of the system followed by an inves- 
tigation of the dissipativc dynamics of the system. 

III. UNITARY EVOLUTION 
A. Schrodinger picture 



and the collective spin of the ultracold gas. The ratio 
between these interaction terms is 

H Zeeman ^ ^ d\fN ^ ^g-j 
^ direct ^O-Nrnag 

However, even if this direct interaction is included, the 
quahtative behavior of the system does not change. We 
must stress, though, that the role of the BEC is essential 
for generating a time independent statistical mixture of 
(dark) entangled states. In the absence of atoms a pair 
of nanocantilevers interacting via the dipole-dipole inter- 
action may well exhibit time dependent entanglement. It 
would, however, be difficult to use or capture these states 
because one would have to, for instance, be able to con- 
trol the direct interaction between the cantilevers such 
that it could be switched off at a chosen time. 

If we assume j22| that the cloud of ultracold atoms can 
be represented by a single quantized mode of frequency 
oja, which can be tuned to resonate with the frequencies 
Wo of the quantized nano-rcsonators, then the dynamics 
of the two nanocantilevers interacting via the ultracold 
gas is well described by the Jaynes-Cummings Hamilto- 
nian [H. 

With the dipole and rotating wave approximations, 
and taking Ka = Kb — k, the Hamiltonian for our sys- 



We shall now consider a unitary evolution of the sys- 
tem of two nanocantilevers interacting via an ultracold 
gas of atoms. As mentioned before, we assume that the 
nanocantilevers have been cooled near to their quantum 
ground states and that only the few lowest-energy Fock 
states of the quantized oscillators are populated. In the 
Schrodinger picture, we restrict the excitation subspace 
of the two nanocantilevers to one, two and three quanta 
of vibrational excitations. 

To study the dynamics in the Schrodinger picure we 
shall assume that the ultracold gas comprises N two-level 
atoms, and that the maximum number of excitations in 
the gas is restricted to one. This is a reasonable assump- 
tion since we have found that in an n-excitation subspace 
with n <ti N, inclusion of higher excitations in the gas 
changes the Rabi frequency without interfering with the 
qualitative behavior of the system. With this assump- 
tion, there are only two global states of the ultracold 
BEC, namely | — J, J) and | — J + 1, J). The action of 
J+ and J_ on these states is 

J+\-J,J) = Vn\-j+i,j) (9) 
J-\-J + l,J) = \/iV|-J,J). (10) 

In the discussion to follow, we denote the collective 
ground state | — J, J) of the ultracold gas by \g) and 
its first excited state \ — J + 1, J) by |e). 
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To a first approximation, we neglect any dissipation 
present in the system and use the Schrodinger equation 
to study the unitary dynamics of the system. In the inter- 
action picture, a general initial state of the cantilever-gas 
system in the one-excitation manifold can be written as 

|*i(t)> = e-^«°*/''[C,,o,iWl5,0,l) (11) 
+ Cg,,,o{t)\g,l,0)+C\o,o{t)\e,0,0)], 

where \g, I, m) (|e, I, m)) represents a state with the ultra- 
cold gas in its ground (excited) state, cantilever a with 
I excitations and cantilever b with m excitations. Using 
([TT|). the Schrodinger equation becomes 



where 



^Ci(t) = /3Ci(0, 



0k 
0k 

K K 



(12) 



and 



Ci(t) 



Ce,0,o(i) 



The coupling matrix /3 has the eigenvalues Aq = 0, 
A±=±i-\/2K, with the corresponding normalized eigen- 
vectors given by 



|Ao) = 



1 

-1 



IA+) = 



|A-> 

1 
1 

72 



-1 
-1 
V2 



(13) 



The eigenstate |Ao) has no excitation in the atomic gas. 
Such a state is termed a dark state or trapped state, 
because its population does not decay even in the pres- 
ence of dissipation in the atomic gas. Dark states have 
been extensively studied in the context of decoherence- 
free subspaccs There have been many proposals 

to prepare qubits in such states and hence protect them 
from decoherence [1^. The population in states orthog- 
onal to dark states undergo dissipation and ultimately 
relax to the ground state, while the population of the 
dark states remains intact. 

If the two nanocantilevers couple to the ultracold gas 
with different coupling strengths then the dark 

state in the one-excitation manifold takes the form 



|Ao) 



1 



(Ka|5,0,l)-K,|ff,l,0)), (14) 



which is a maximally entangled state of the two nanocan- 
tilevers for Ka = Kf,. It is interesting to see from Eq. (1141) . 



that by tuning the interaction strengths Ka and Kb, adi- 
abatic state transfer can be achieved between the 
two cantilevers. This can be realized, for instance, by 
instead using electromagnets, making it possible to vary 
the coupling strengths between the vibrational mode of 
each nanocantilevcr and the spin of the ultracold Bose 
gas. The fabrication of such electromagnets on nanocan- 
tilevers is obviously a challenging task, but recent ad- 
vancements in nanofabrication techniques suggest that it 
might be possible in the future [11]. We will further dis- 
cuss the relevance of dark states for our scheme in the 
next section, where we will study the dissipative dynam- 
ics of the system. 

Equation PT|) can easily be solved, giving the time- 
evolved wave function in the onc-excitation manifold. 
With Cg^ifi{0) = 1, the system evolves to 



1 



cos(\/2Kt) 
cos(\/2Kt) — 1 



(15) 



The probability for each of the three basis states to be 
occupied is plotted as a function of time in fig.[H The ex- 
citation is reversibly transferred between the cantilevers 
and the ultracold atomic gas. We shall explicitly quan- 
tify the entanglement present in different excitation man- 
ifolds in section |Vl but note already here that the states 
|(7,1,0) and 0,1) exhibit time-varying entanglement. 
It is also worth noting that the amount of entanglement 
between the two nanocantilevers strongly depends on the 
initial state of the system. The amount of entanglement 
generated between the two nanocantilevers is enhanced 
if the initial excitation lies in the ultracold gas. 

To study the dynamics when more than one excita- 
tion arc present in the system, equation (jlip can readily 
be generalized. In the interaction picture, general state 
vectors for the system in the two- and three-excitation 
manifolds take the form 



l*2(i)) = 



H„t 



l*3(i)> 



J=0 



■^Ce,j,l-jit)\e,j,l- j) 



E^9,..3-,(t)|5,J,3-j) 



(16) 



(17) 



Using pB)) and pT)) , the Schrodinger equation reduces to 
sets of coupled differential equations which can be easily 
solved. A typical solution for a three-excitation subspace 
is shown in fig. [31 where occupation probabilities are 
plotted as functions of time for the different basis states. 
The result is a time-varying entanglement between the 



5 




FIG. 2: (Color online) Time evolution of the occupation prob- 
ability for (a) the states Ig, 1,0) (dashed), \g,0,l) (dotted), 
and |e,0, 0) (solid) with the initial condition Cg,i,o(0) — 1, 
and (b) for the states Ig, 1,0) and \g,0,l} (dashed, identical) 
and |e, 0, 0) (solid) with the initial condition Ce,o,o(0) — 1. 



FIG. 3: (Color online) Occupation probability, as a function 
of time, for the states (a) \g,3,0) (dashed), \g,2,l} (thick 
dashed), \g,l,2) (dotted), and |g,0, 3) (solid); (b) |e,2,0) 
(dashed), |e, 1, 1) (dotted), and |e, 0, 2) (solid), with the initial 
condition Cg,3,o(0) = 1. 



states in each excitation manifold. It turns out that the 
system of two nanocantilevers exhibits entanglement in 
higher excitation subspaces, even for mixed initial states, 
such as a thermal state. Thus thermal effects do not nec- 
essarily degrade the entanglement of a quantum system. 



B. Heisenberg picture 

The Schrodinger and Heisenberg pictures are two 
equivalent ways of studying the evolution of a quantum 
system. After studying the system in the Schrodinger 
picture, we now present a compact way of studying the 
unitary evolution of the system in the Heisenberg picture. 
This method, apart from giving more insight into the dy- 
namics of the system, also allows us to exactly solve the 
dynamics of the system for arbitrarily many excitations. 

Using the Dicke model introduced in sectionlHl we have 
— J < A/ < J. Since most of the atoms in the ultracold 
condensed in the ground state, and is large, we 
can approximate [J+, J_]=2Jj « —N. An equivalent way 
of invoking this approximation is through the Holstein- 
Primakoff (HP) transformation This is essentially 



a mapping from angular momentum spin operators to 
bosonic operators. Under the HP transformation, the 
collective angular momentum spin operators transform 
as 

J+ = cWN-cU (18) 
J_ = c^N ~ etc (19) 
J, = ctc-7V/2, (20) 

where c and & are the bosonic lowering and raising op- 
erators which satisfy [c, c^] = 1. In the large N limit, 
Eqs. nil)- (HOI) reduce to J+ « ct//V, J_ « c//V and 
« etc — A^/2. Under this approximation the Hamilto- 
nian ([5]) becomes 

— =a;o(c^c-A^/2 + ata + 6tS) + [K(a + 6)ct + /i.c.]. (21) 

We now introduce pj-(O) and g^(0) as two sets of 
bosonic creation operators for collective excitations. By 
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choosing 



qHO) 



i[at(0) + 6t(0)]±-Lct(0) 
1 
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[at(0)-6t(0)], 



the Hamiltonian (|2T|) becomes diagonal, 



(22) 
(23) 

(24) 



where Q± = wq ± v^k. The time-evolved creation oper- 
ators are 



qt(t) = e'"°*gt(o). 



From (pSj) and (^5)) we obtain 

at(i) = i[p+tW+P-^W + V2g^W] 
= i[pVW+P-^(i)-V2gt(i)] 



(25) 
(26) 



(27) 
(28) 



These time evolved creation operators (t) and 6^ (t) can 
conveniently be used to describe the unitary time evolu- 
tion of the system of two nano-cantilevers and an ultra- 
cold gas for arbitrary initial conditions. Evolution of a 
state with cantilever a initially in its nth excited state 
can easily be determined from |vl'„(i)) = {d^t))"\g,0,0) . 
For instance, if n = 1, we obtain 



\^,{t))=dHt)\g,0,0) = l 



1 + cos{V2Kt) 
cos{V2Kt) - 1 



(29) 



where the basis vectors are the same as for Eq. ([T5t . 
obtained in the Schrodinger picture. The results are 
identical apart from a phase shift on the last compo- 
nent; this is a general feature of the Schrodinger vs. the 
Heisenberg picture. The dark states of the system are 
seen to correspond to excitations of the mode labeled 

so that a dark state with n excitations is given by 
\^q^ri{t)) = (gt(t))"|g,0,0). This way, dark states corre- 
sponding to different numbers of excitations in the sys- 
tem can easily be computed. 

Before concluding this section, we shall briefly point 
out the similarities and differences in the Schrodinger and 
Heisenberg approaches. To simplify calculations without 
compromising physical insight, in the Schrodinger picture 
we restricted the maximum number of excitations in the 
atomic gas to one. In the Heisenberg picture, this con- 
straint is relaxed, and the system of two nanocantilevers 
interacting with the ultracold BEC is effectively treated 
as three coupled harmonic oscillators. In reality, however, 
the number of excitations in the atomic gas is limited by 
the total number of atoms N. The Heisenberg solution 
is therefore valid as long as the total number of excita- 
tions is less than the total number of atoms N, so that 



the gas is effectively equivalent to a harmonic oscillator. 
The Schrodinger and Heisenberg approaches give identi- 
cal results if this holds, and if more than one excitation 
is allowed in the gas in the Scrodinger picture. 



IV. DISSIPATIVE DYNAMICS 

So far we have neglected decay channels in our system, 
but with any real physical system, dissipation is of course 
inevitable. Dissipation can occur in the atomic gas with a 
decay rate F. or by thermal decay of the nanocantilevers 
with a decay rate 7. For a nanocantilever with resonant 
frequency ujo/2tt = 1 MHz and quality factor Q = 10^ 
[soj . one obtains 7 ^ 1 Hz 0], whereas F is largely gov- 
erned by spin flips due to collisions or stray currents in 
the magnet or the cantilever [3l|, [s^. We have chosen 
F = 10 Hz in the examples below, which allows us to 
neglect any direct cantilever dissipation and consider the 
much more rapid decay of the ultracold gas as the only 
dissipation channel. We note also that if the ultracold 
gas loses an excitation, then it is very unlikely that a 
cantilever can pick it up, owing to their very small solid 
angles as seen from the ultracold gas. 

Under the Born-Markov approximation, evolution in- 
cluding dissipation [s^ is well described by a Lindblad- 
type master equation of the form 



where p is the density matrix of the system, and 



(30) 



Cp = -{2J^pJ+ - J+J-p - pJ+J-) (31) 

is the Lindblad operator. With at most one excitation 
present in the system, Eq. ([3T|) can be solved by evalu- 
ating 



[Cp) 



■1,3 



(32) 



, fc=0.i=0 



where |0) = |5,0,0), |1) = |e,0,0), |2) = |5,1,0) and 
|3) = \g,Q, 1). In this basis, Eq. ([50]) transforms to a set 
of coupled differential equations. 



^ p+Uh,p\ = 



dt 



NT 



-2pi,i po,i 

Pi,o 2/91,1 Pi. 2 

P2,l 

P3,l 





Pi. 3 





(33) 



A numerical solution of Eq. ([33]), for the initial condition 
(^9,1,0(0) = 1, is shown in Fig. |4l In the steady state, 
the system relaxes to a statistical mixture of the ground 
state and the dark state with equal probability. Hence 
dissipation-assisted time evolution leads to a long-lived 
entangled state of two nanocantilevers. 

Following a similar series of steps as outlined above, 
we have solved for the dissipative dynamics with two 
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FIG. 4: (Color online) Time evolution of the occupation 
probability for the states \g,0,0) (solid), |g, 1,0) (dotted), 
\g,0,l) (thick dashed) and |e,0,0) (thin dashed). Initially 
Q,i.o(0) = 1. 



and three excitations initially present in the system, and 
have found similar behavior. For an initial state where 
(^9,2,0(0) = 1, a numerical solution of Eq. pop is shown 
in fig. [5l The system relaxes to a statistical mixture of 
dark states with zero, one, and two excitations. With 
the chosen initial state, the probabilities for these dark 
states to be occupied are 0.25, 0.5 and 0.25 respectively. 
For three excitations in the initial state (not shown here), 
the steady state is a mixture of dark states with zero, 
one, two and three excitations, with respective statisti- 
cal weights of 0.125,0.375,0.375 and 0.125 for an initial 
state Cg^3fi{0) = 1. 

The robustness of the scheme lies in the fact that en- 
tangled states of nanomechanical systems may be gener- 
ated even if these are initially found in thermal states. 
This will be further investigated in future work. Here, 
we study a simple example. We solve for the dynamics 
of the system with one of the cantilevers and the ultra- 
cold gas prepared in their respective ground states. The 
other cantilever is in a mixture of zero, one, two and three 
excitations, with a distribution that approximates that of 
a thermal state. The solution of the master equation is 
shown in fig. [51 and the resulting steady state is entan- 
gled. To summarize this section, we have shown that 
two nanocantilevers, interacting with the same dissipa- 
tive cloud of ultracold atoms, will become entangled with 
each other. The entangled state is immune to dissipation 
in the ultracold gas, and hence opens the possibility of 
achieving long-lived entangled states of nanomechanical 
systems. 



V. ENTANGLEMENT MEASURES 

After showing that it is possible to generate long- 
lived entangled states of two nanocantilevers, we will 
now quantify the entanglement between these. In sec- 
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FIG. 5: (Color online) Occupation probabilities, as a function 
of time, for the states (a) \g,0,0) (solid), |g,0, 1) (dotted), and 
|g, 1,0) (dashed); (b) |g,0,2) (solid), |g, 1, 1) (dotted), and 
\g,2,0) (dashed). Initially Cg, 0,2(0) — 1. Excitations in the 
ultracold gas decay so quickly that the probability for states 
containing such excitations to be occupied are much smaller 
than the probabilities shown here. 



tion IIII Al we used the Schrodinger equation to study 
the unitary evolution of the system of two cantilevers in- 
teracting with an ultracold BEC with one, two and three 
excitations present in the system. The time evolved state 
of the two nanocantilevers is a mixed state found by trac- 
ing over the atomic gas. 

To quantify the entanglement in a bipartite system in a 
mixed state, we use the Peres criterion [3j|. We compute 
the negativity for the reduced density matrix, defined as 



max{0, — Xi) 



(34) 



where Xi is the sum of all the negative eigenvalues 
of the partially transposed reduced density matrix [ssj . 
The negativity Af is not only easy to compute but has 
an added advantage of being an entanglement monotone 

Using the time evolved wave functions obtained in 
IIII Al the negativity A/", as defined in Eq. ([34)) . is com- 
puted and is shown in fig. [3 For a bipartite system of two 
qubits, Af lies between and 0.5, with A/" = for separa- 
ble states and Af = 0.5 for maximally entangled states. 
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FIG. 6: (Color online) Occupation probabilities, as a function 
of time, for the states (a) \g,0,0) (dashed), (b) Ig, 1,0) (dot- 
ted), |g,0, 1) (solid) , (c) \g,2,0) (dashed), |g, 1, 1) (dotted), 
|p,0,2) (solid) , (d) |5,3,0) (solid), |g,2, 1) (thick dashed), 
\g,l,2) (dotted) and \g,0,3) (thin dashed). Cantilever 'a' is 
initially in a mixed state of zero, one, two and three excita- 
tions, approximating a thermal state with {rithermai) = 0.3. 
Excitations in the ultracold gas decay so quickly that the 
probability for states containing such excitations to be occu- 
pied are much smaller than for the other states considered 
here. 



For a pure state of two qudits, with k terms in its Schmidt 
decomposition, Q < N < [k — l)/2. As shown in fig. [3 in 
the one-excitation subspace the two nanocantilcvers are 
entangled at all times, except at certain instants when the 
negativity becomes equal to zero. Interestingly, in higher 
excitation subspaces, the two nanocantilevers are always 
entangled for t > Q. In fig. [7] we have also compared N 
for two different initial conditions in the one-excitation 
subspace. As is clear from fig. [7l if the initial excitation 
lies in the ultracold gas, then maximally entangled states 
of the two nanocantilevers occur at certain times. One of 
the main features of our proposed scheme is that entan- 
gled states of the two nanocantilevers are generated even 
for an initial classical state. As an illustration, we plot 
in fig. |8]the negativity as a measure of entanglement for 
an initial separable mixture of the first three excitation 
subspaces. A non zero value of negativity in fig. [S] again 
points to a finite entanglement being present between the 
two nanocantilevers. 

As discussed in the previous section, in the presence of 
dissipation, a system of two nanocantilevers coupled to 
an ultracold gas of atoms relaxes to a statistical mixture 
of various dark states. The dark states are pure entangled 
states of the two nanocantilevers. For pure states, it is 
possible to quantify entanglement using the participation 
ratio C 37[, defined as 



1 



(35) 



where Pi are the coefficients in the Schmidt decomposi- 
tion of the pure state \^) of the bipartite system. The 




FIG. 7: (Color online) Degree of entanglement, as measured 
by the negativity defined in Eq. (|34[) . for a system of two 
nanocantilevers interacting with a dissipation free ultracold 
gas, in the one-ecxitation (thick solid), two- excitation (dot- 
ted) and three-excitation (dashed) subspaces, with all the ex- 
citations initially present in one of the cantilevers. Also, for 
comparison, the negativity is presented for the case when the 
initial excitation is in the gas for the one-excitation subspace 
(thin solid). 
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FIG. 8: (Color online) Degree of entanglement, as measured 
by the negativity for a system of two nanocantilevers interact- 
ing with a dissipation free ultracold gas for an initial mixed 
state of the first three excitation subspaces with average ther- 
mal occupancy of 0.3. 



quantity ^ lies between 1 for separable states and k 
for maximally entangled states. For the dark states in 
the one-, two- and three-excitation subspaces we have 
Ci,2,3 = 2, 2.67 and 3.2 respectively. Although Cs > C2 > 
Ci, only the dark state in the one-excitation subspace is 
a maximally entangled state of the two nanocantilevers. 
Nonetheless, in higher excitation subspaces there is a 
finite degree of entanglement present between the two 
nanocantilevers. 



VI. CONCLUSIONS 

We have proposed a scheme to generate entanglement 
between two nanomechanical systems through indirect 
interactions mediated by an ultra cold gas of atoms. The 
excitation energy of the two-level atoms is assumed to 
have been tuned to coincide with the energy quantum of 
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the fundamental flexural mode of the two identical can- 
tilevers. Hence single excitations have the same energy 
in the atomic system as in the cantilevers and we have 
studied the dynamics of systems with a total of one, two 
and three excitations present. In the presence of dissi- 
pation in the atomic gas, the two nanocantilevers relax 
into a statistical mixture of various dark states of the 
system. This opens the possibility to achieve long-lived 
entangled states of nanomechanical systems which are 
only limited by the dephasing time of the cantilevers. 
The atomic decay channel is important for the creation 
of time independent entangled states, provided that the 
decay rate of the atoms is low enough to enable a strongly 
coupled regime and high enough so that the steady state 
is reached before the coherence of the cantilevers decays 
away. 

With respect to the preparation of an initial state it 
has been shown in [ssf that if a two level system is cou- 
pled to a classical field and a quantum field, prepared in 
its vacuum state, and if the two coupling strengths can be 
tuned independently of each other then atom-field inter- 
actions can force the vacuum state of the field to an arbi- 
trary superposition of the Fock states. This can be used 
for the initial state pre par ation of the nanocantilevers. 
Also, in a recent work [39[ a possibility of transferring 
also non-gaussian states from an optical field to a me- 
chanical oscillator has been outlined. 

A challenging aspect of any scheme involving entan- 
glement generation between "macroscopic" mechanical 
systems is the actual experimental detection of entan- 



glement. A promising scheme suggested in [40| employs 
a quantum non demolition (QND) measurement of the 
excitation of a mechanical oscillator by using the an- 
harmonic coupling between two beam bending modes, 
which uniquely shifts the resonant frequency of the read- 
out oscillator. Another recent proposal is to determine 
the Wigner function of a mechanical resonator by cou- 
pling the vibrational mode of the oscillator to a three 
level "detector" atom, which in turn is coupled to a pair 
of optical fields, inducing Raman transitions between its 
ground and excited state [4lj. It has been shown that 
the probability for the "detector" atom to be found in the 
ground state is a direct measure of the Wigner character- 
istic function of the nanomechanical oscillator as long as 
dissipation in the system can be kept under control. En- 
tanglement detection for nanomechanical systems is thus 
experimentally challenging, but may be possible using 
techniques such as these. 
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